Interference of parametrically driven one-dimensional ultracold gases 



Susanne Pielawa 
Physics Department, Harvard University, Cambridge, MA 02138, USA 
(Dated: February 3, 2011) 

We theoretically analyze interference patterns of parametrically driven one-dimensional ultracold 
atomic gases. By modulating the interaction strength periodically in time, we propose to excite 
collective modes in a pair of independent one-dimensional gases at energies corresponding to the 
drive frequency. The excited collective modes lead to spatial oscillations in the correlations of the 
interference pattern, which can be analyzed to obtain the sound velocity of the collective modes. 
We discuss both bosonic and fermionic systems, and how such experiments could be used to probe 
spin charge separation. 



I. INTRODUCTION 

One-dimensional systems show many interesting prop- 
erties like charge fractionalization and separation of 
charge excitations from spin excitations. There is exper- 
imental evidence for these properties in various systems 
[THZ]. Recent progress in trapping and cooling of cold 
atoms, has made it possible to study one-dimensional sys- 
tems and their properties in both fermions and bosons j8j- 
113] . One of the primary tools in measuring properties of 
cold atoms are interference experiments. 

Surprisingly, it has been shown that interfering two in- 
dependently created Bose Einstein condensates (e.g. by 
releasing the confining potential and letting them expand 
and overlap) gives rise to spatially periodic patterns in 
each experimental shot. As there is no well defined phase 
between the two condensates, the maxima and minima 
of the pattern are at different positions in each shot, so 
that the average over many shots, corresponding to a 
quantum-mechanical ensemble average (p(r)) = 0, con- 
tains no interference fringes. The key to this effect is that 
we can decompose the state of the two condensates, each 
containing a definite number of particles, into a superpo- 
sition of states of well defined relative phases. Measur- 
ing the interference pattern collapses the wavefunction 
of the combined system to a state with definite phase, 
and thus produces a periodic pattern. To theoretically 
study this effect, we need to study the density-density 
correlation function (p(r 1 )/5(r 2 )) instead of the average 
density, since the quantum mechanical average of the 
density formed in interference experiments samples all 
possible relative phases and thus hides the interference 
pattern. The above arguments hold true even for an in- 
terfering pair of one-dimensional systems, with one im- 
portant modification: as there is no long range order, 
the maxima formed in the interference pattern are wavy 
lines instead of straight lines. This waviness contains in- 
formation about the correlation function within each Id 
system [THUS]. 

In this paper, we propose an experimental approach to 
create and directly probe excitations in one-dimensional 
ultracold atom systems of both bosons and fermions. 
In particular, for systems composed of mixtures of two 
species, our method allows us to see separation to spin 



and charge modes. Our approach is to start with a pair 
of one-dimensional systems, parametrically drive excita- 
tions by temporally periodic modulation of the parti- 
cle particle interaction strength, and finally analyze the 
interference pattern. Experiments involving parametric 
driving to create excitations in one-dimensional fermionic 
gases have recently been proposed in Refs [THJ [T7] . In 
this work, in addition to fermionic gases, we study also 
bosonic gases, which are of experimental interest. 

In interference experiments with one-dimensional sys- 
tems of cold atoms, the gases are released from their traps 
and allowed to expand. We label the longitudinal direc- 
tion as the z the axis, and the two systems lie in the x-z 
plane, see Fig. [T] Due to the strong transverse confine- 
ment, when released from their traps, the atoms expand 
mainly in the transverse directions. Therefore, we shall 
neglect the motion of atoms in the longitudinal direction 
during the expansion. Once the gases have expanded to 
a size much larger than the initial distance d between the 
two systems, an absorption image is taken. We neglect 
the slowly varying envelope and the expectation value of 
the density in case of independently created quasicon- 
densates is a constant, (p(x,z)) — ptofPo, where po is 
the line density of each quasicondensate. The density- 
density correlations function, however, contains an oscil- 
lating term at wave vector Q = m is the particle 
mass, t the time of expansion: 



(p(n)#ra))oc(p8± cos(AxQ)J, (1) 

where rj stands for (xi,Zi), Ax = x\ — xi, the upper 
(lower) sign applies to bosons (fermions). ip(z) denotes 
the particle annihilation operator in either of the equiv- 
alent initial Id systems. To derive Eq. ([I]) we have ne- 
glected the motion of atoms along the longitudinal direc- 
tion during expansion; assumed that the two systems are 
initially identical and independent; assumed large expan- 
sion times (e.g. the size of the overlapping clouds is much 
larger than their initial separation d). See appendix A 
Ref [T5] contains a more detailed discussion. 

Typically, experimental data is analyzed by looking 
at the the integrated density-density correlation func- 
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We define 



quasicondensate 1 



quasicondensate 2 



p(x,q) = / e tqz p(x, z)dz 



(2) 



and then consider the quantity (p(x, q)p(0, —q)) which is 
given by 

(p(x,q)p(0,—q)) oc cos(xQ)L / dzcos(qz) (^(z)V'(O)) 
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Thus, by Fourier transforming the interference pattern, 
one can obtain the Fourier transform of the correlation 
function of a single ID gas, 



\A(q,T)\ 2 := L / dzcos(qz) (ft(z)i>(0)) 



(3) 



The quantity ^(g,! 1 )! 2 is also a function of the total 
driving time T (we have suppressed this time argument 
T in ip and p in the above equations). This analysis can 
in principle be done for fermions and bosons. 



FIG. 1: (Color online) Proposed setup: two independently 
created one-dimensional systems of ultracold atoms are first 
parametrically driven by changing the inter-atomic interac- 
tion strength periodically in time (e.g. by changing the trans- 
verse trapping frequency), at time T the driving stops and the 
systems are released from their traps and allowed to overlap. 
An absorption image is taken, showing wiggely fringes along 
the z axis. These fringes contain information about the two- 
point correlation function within each gas (before expansion) 
and can be analyzed to obtain the sound velocity of the col- 
lective modes. 



tion HUES], 



dz 1 dz 2 (p(x, Zi)/3(0, z 2 )) oc p L + (\A L \ ) cos(Qx) 



where L is the length of 
interference amplitude Al 

L J L dz 



and 
(\Al\ 2 ) 



the 



the integration, 
is given by 

One can now define the inter- 
\Al\/pqL, which decays as the 



(VW(o)> . 

ference contrast Cl = 
integration length L is increased, encoding the decay of 
the correlation functions in the one-dimensional systems 

EH- 

Our approach is to study the interference pattern 
of one-dimensional gases that have been parametrically 
driven out of equilibrium. The driving is done by a tem- 
porally periodic change of interaction strength for a cer- 
tain time T before the gases are released from their traps. 
To extract information from the interference pattern, we 
suggest to take the Fourier transform along the z-axis, 
and then study correlations for that Fourier transform. 



Our main results are: (1) We find that parametric drive 
at frequency fl creates pairs of excitations with total mo- 
mentum zero. Therefore, the quantity shows a 
resonance peak at the wave vector corresponding to half 
the driving frequency, q = H/2ui (and, in principle, as- 
sociated multiples) , where m is the sound velocity of the 
corresponding mode. Such an analysis can be used to find 
the sound velocity of the ID bose system. (2) For a two 
component Bose system, we find two primary peaks cor- 
responding to "spin" and "charge" velocity. (3) Similarly 
one could analyze an interference pattern of a driven ID 
system of two component fermions, where one expects 
two peaks corresponding to spin and charge excitations. 
For fermions, however, there are some complications due 
to the rapid decay of the correlation functions, which we 
will address in Section llVl 



The paper is organized as follows: In Section [II] we in- 
troduce the formalism and describe the relation between 
a driven Id gas and a driven simple harmonic oscilla- 
tor. In Section ITO1 we calculate the correlation functions 
of driven systems and show how information about the 
driving can be extracted from the interference pattern. 
In Section |IV| we address the problem of observing spin 
charge separation for fermionic Id systems. 
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II. PARAMETRIC DRIVING AND LUTTINGER 
MODE SQUEEZING 

A. Description of a Id system: Luttinger liquid 

One-dimensional systems of interacting bosons or 
fermions can in general be described by Hamiltonian 



a condition that can be realized in experiments using hy- 
perfine states \F = l,mp = —1) and \F = 2, mp = +1) 
of 87 Rb [H[5TJ[52]. We label 9± thc interaction strength 
between different species. In this case the Hamiltonian 
separates into a "charge" part with g c = g|| + g± and 
"spin" part with g s = g\\ — g±. The ground state corre- 
lation function becomes 



H = 



L r L 

g aa \z - z' 



(4) 







z 



(I 



where g{z) is the interaction strength, and i/>'{z) creates 
a particle at point z. Low energy excitations may in 
turn be described by the Luttinger liquid Hamiltonian 
which in momentum space is given by a sum of harmonic 
oscillators 



H = 



2^ 

q^O 



uK 



where n and <p are conjugate fields: 



<f<Pq0-q 



iS a 



(5) 



-q 1 1 



K and u are the Luttinger parameters, u corresponds to 
the sound velocity, and K is a dimensionless parameter. 
For a finite system of length C we quantize the momen- 
tum g = ^-n. We shall first consider the case of spin- 
less bosons and then generalize to a two component Bose 
mixture and fermions with spin. For a weakly interact- 
ing one component Bose system (K > 1) the Luttinger 
parameters are given by 



K = why/po/gm, 



(6) 



where po is the ID line density, and g is the interaction 
parameter, related to the three dimensional scattering 
length a s through g — 2hi>T a >s [HE 120]) where is the 
transverse trapping frequency. The relation between the 
physical boson creation operator and the Luttinger fields 
is given by 



V> T (z) = \Jpo+Tl(z)e 



-iip(z) 



(7) 



and the correlation function in the ground state has a 
power law decay 



(f(^(0))o«ft 



(8) 



where ^ = 2K/ po is a short distance cutoff. Let us now 
assume that there are two species of bosons (e.g. two 
hyperfine states or two different atoms), labelled j" and 
4-. We further assume that the interaction between two 
atoms of the species labelled up and the one between two 
atoms of the species down are the same, g^ = g^ = , 



(V t (z)^(0)) O (xpo ( 



(9) 



where K c (K s ) is related to g c (g s ) through Eq. 



Next let us consider a fermionic system of two compo- 
nents labelled t an d 4-- I n t ne case of ultracold atomic 
gases, two fermions of the same species cannot interact 
via s-wave scattering (g^ =0). In thc following we as- 
sume (^-interactions and denote by g the interaction con- 
stant between two fermions of opposite spin. The Hamil- 
tonian separates in a charge part H p and a spin part H a . 
Both charge and spin part have the form of Eq. ([5| with 
thc Luttinger parameters given by |23j 



u v K v = vf, 

".- ( - ■ 9 \ 



VF 1± 



TTVp J 



(10a) 
(10b) 



here v = p, a and the upper (lower) sign applies to the 
charge (spin) mode. The single particle annihilation op- 
erators are given by ^{z) = ipf l u(z)+ipf t L(z), where the 
index R (L) stand for a right (left) moving particle, and 



Poe - lkFZ e 



p Qe ik F z e - i ( t Pr( z )+ t( z )) 



i(-<p f (z)+0f(z)) ^ (11a) 
(lib) 



analogously for 4-. The fields which decouple the Hamil- 
tonians in H p and H a are given by 



1 

71 
i 

71 



(12) 
(13) 



and same for 9 U , which enters in the Luttinger liquid 
Hamiltonian as W6 u (z) — ■kT\ v (z). 

The spin part of the Hamiltionan has an additional 
sine Gordon term, arising from backscattering between 
fermions with opposite spin (<7i_l term). However, this 
term renormalizes to zero |23j if K a > 1 (repulsive in- 
teractions), otherwise it creates a gap in the spectrum of 
the spin Hamiltonan. In this paper we concentrate on 
repulsive interactions and the sine-Gordon term is not 
relevant. 

The single particle correlation function in the ground 



4 



state 



differential equation 



(-(/>t 0)^(0)) = 2 Po cos(k F z) 



6, 



(14) 



is oscillating at the length scale of particle spacing. 



B. Parametric driving and squeezing 

Experimentally the interaction parameter g can be 
changed by either changing the scattering length a s , or by 
changing the transverse confinement. Let us now imagine 
a Id system that is parametrically driven by changing g 
periodically in time for a certain time T, 



g(t)=g(l-S'oos(Qt)), 



(15) 



where g is the average value of interaction and Q is 
the driving frequency. Since the Luttinger Hamilto- 
nian Eq. ^ is just a collection of harmonic oscillators, 
this situation can be mapped to a collection of para- 
metrically driven simple harmonic oscillators via P q = 
— y/uKq 2 /TTifg and Q q — ^J-n juKq 2 i\ q . We note that 
this mapping is time-independent for both the bosonic 
and the fermionic case, as uK remains constant. Each 
pair of modes (q, —q) is then described by 



H q {t) = ^P- 



[l-8cos(ttt)]Q q Q- 



(16) 



with 



and u 2 



9Po, 



uq, 



5' for bosons and u 2 



(17) 



,(1 ± -2- 



±6'- 



-r= for fermions. It is well know that a para- 

metrically driven quantum harmonic oscillator can be 
mapped to a static reference system 24J. If initially in 
the ground state, the final state after a driving time T 
is, up to an overall phase factor, given by S (£ q (T)) |0) 
where 



Sfo(T)) = exp e q (T)a q a_ q - U T H^- g 



(18) 



is the squeezing operator, and a q (a q ) is the annihilation 
(creation) operator of one excitation in the static refer- 
ence system, 



Wref(<7) 



Qa 



p 



(19) 



We choose our reference system to have the eigenfre- 
quency of the system at t = 0, u> re f(q) = ui q yl — 5. The 
squeezing parameter £ 9 (T) can be expressed in terms of 
solutions of the classical parametrically driven harmonic 
oscillator, namely solutions to the following second order 



z"(r) + 



[1 -Jcos(2t)]z(t) =0, 



(20) 



which we have expressed in dimensionless variables, and 



We define r\ = ^r, and 



??rof 



r]y/l — S. Let 



X (r) and y(r) be solutions of this equation with X(0) 
y'(0) = 1 and X'(0) = y(0) = 0, and define z{t) 
X(t) + irj rc fy(r), then ^ q is defined through 



4cosh 2 (|£ 9 |) = |z| 2 + |zfA7 r 2 of + 2 



(21) 



and an expression for its phase, £ = |£| e l1? , can be found 
through = d u + § v , where \u\ = cosh |^|, and [24] 



\u\e^ 
\v\e^ 



2 



z(t) 
z(t) + 



iz{t) 
iz(t) 



(22a) 
(22b) 



Physically relevant quantities can be expressed in terms 
of the above, e.g. average number of excitations (a' a) = 
sinh 2 |£|, and uncertainties in position and momentum 
variables: 7y rcf ^Ai 2 = \z\ /2 and r] le f^Ap 2 = \z'\ /2. 



We note that Eq. ( |20| is Mathieu's equation and its 
odd and even solutions X(r) and y(r) are referred to 
as Mathieu S and Matheiu C functions, because of their 
resemblance to sinusoidal functions in certain parameter 
regimes. Depending on parameters rj and S, this equation 
has stable, i.e. (quasi-)periodic solutions or unstable i.e. 
exponentially growing and decaying solutions, see Fig. [2] 
On resonance 77 ~ 1 the squeezing parameter grows as 

|£|«£rn. 

To summarize the above, the final state of a paramet- 
rically driven Id system is fully described by the set of 
squeezing parameters £ g , which in turn is given by solu- 
tions of Eq. ((201). 



III. 



DETECTING EXCITATIONS VIA 
INTERFERENCE 



We now show how the excitations of the driven con- 
densate change the correlation function and the interfer- 
ence pattern. While the parameter K can be extracted 
from the ground state interference pattern, interference 
of driven systems could be used to measure the sound 
velocity u. Given a set of squeezing parameters £ g we cal- 
culate the correlation function for bosons and fermions 
and obtain {\A q \ ), a quantity that can be extracted from 
the interference pattern in experiments. 
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b) 





n =i.2 
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FIG. 2: (Color online) a) Stability diagram of Eq. (20 1, as a function of momentum q (in units of ||) and relative driving 
amplitude 8. Gray shaded regions are resonant, where the number of excitations in that mode grows exponentially in time. 
White regions are off-resonant and the number of excitations is a (quasi-) periodic function in time and thus bounded from 
above, b) Absolute value of squeezing parameter |£| as a function of time for different values of r\ = -S 2 -, at fixed S = .1. In 
the unstable region, £ grows linearly in time. The first resonance occurs at 2uq — f2, when the driving frequency is twice the 
frequency of an excitation in that mode, since excitations are created in pairs (— q, q). 



a) 



b) 




in gound state 




z (units of 2u/0) 



FIG. 3: (Color online) Driven one-dimensional Bose gas: (a) Observable quantity (\A(q,T)\ 2 ) , divided by its ground state 
value, as a function of wavevector q and number of oscillations complete during driving time T. At T = the system is in 
the ground state. As driving continues a resonance peak corresponding to the drive frequency and sound velocity emerges and 
grows in time. The relative driving amplitude here is 8 = .1. (b) Correlation function after driving time T = 20^. 



A. Bosons: one and two component systems 



quadrature (ip q ip^ q ) t due to the excitations, 



The correlation function can be written in terms of the 
quadratures (<p q <p- q )t as 



0, 



-j, fd( V ) sin 2 (>4)d n 



'(23) 



where z = ~z is the normalized dimensionless position 



variable. The function d{rf) is defined as the change in 



«<f< 



~2iT 



((<Pq<P-q)t ~ (<P q <P-q)o) (24) 



and can be expressed in terms of the squeezing parame- 
ters = r^e 1 "® 11 



1 



(cosh(2|£,|)-cosVinh(2|£,|)-l). 



(> 



We calculate the correlation function and (\A q \ ) of a ID 
bose gas after it has been parametrically driven. The 



integral in the exponent of Eq. ( 23 ) is evaluated numer- 



ically. Fig. [3] shows a result, as expected there is a peak 
in \A q \ 2 at a momentum corresponding to the driving 
frequency. 

There are two modes in a two component bose system 
with different sound velocities. The parametric driving 
excited each of the modes independently. This creates 
beats in the correlation function, which can be seen as 
peaks in A q : there are then two peaks, one correspond- 
ing to each mode, see Fig. |4j This setup can be used 
to probe spin-charge separation in continuous bosonic 
systems. It has previously been proposed to use con- 
tinuous bosonic onc-dimensional systems to probe spin- 
charge separation [TU] . 



B. Two component Fermi systems 

We first calculate the correlation function for the 
driven fermi gas. Note that both fields appear in the 
exponent of single particle operators, see Eq. (11). We 
find for each mode 



= poe' tk - z e-'M (25a) 



I(z) 



, ■ 2 ( Vz 
dri sm [ — 



1 

K 



K+^} smh 2 (r ?) )- 



1 

K 



cos d v sinh(2r 



(25b) 



The above formula is valid for spin-polarized fermions, 
and also for spin and charge part of a two component 
system of fermions. In that case, we just have to replace 
the constants (e.g. interaction strength g and what fol- 
lows from it) by the corresponding value for the spin- or 
charge part. 

As a first approximation, let us neglect the second part 
of I(z) with respect to the first part. There are two rea- 
sons for that: first, for fermions the Luttinger parameter 
K is close to one (in fact, for the spin part it renormal- 
izes to K a = 1), second, only the second term contains 
the phase of the squeezing parameter cos$ 9 , which will 
average out if the systems are left to evolve after driving 
and before interfering. We then have 



I(z) fa const — I K 



1 

K 



drj cos (rjz) - sinh 2 r, ; (26) 



does not give any peaks. In the next section we dis- 
cuss why and how to resolve this issue. Using parametric 
driving and analyzing the momentum distribution of cold 
fermionic gases to study spin-charge separation has been 
proposed in [161 117) . Here we propose to analyze an in- 
terference pattern. 



IV. OBSERVING SPIN CHARGE SEPARATION 
FOR FERMIONS 

In Fig. [4] we see two peaks in |A g | 2 for the two compo- 
nent bose system, corresponding to the sound velocities 
of two modes. In this section we show how one could 
observe spin charge separation for fermions in a similar 
way, however, for fermions there are some complications. 
The fast decay of correlation functions coming from the 
factor (1/K + K) in the power law exponent poses a 
problem when calculating (\A q \ ), as this quantity will 
then be dominated by contributions from close to the 
origin z = 0, preventing the Fourier transform from re- 
solving the beats in the correlation function. We can 
solve this problem by using a simple trick to avoid the 
region around the origin (similarly as used in 25 J. Wc 
need to extract 



cos{zq) (4>\z)iP{0) (27) 



from the interference pattern, with z cnd > z bcg > 0. 

This can be realized by looking at correlations of 
Fourier transformations of density of different regions in 
the interference pattern. For this purpose we define 



p 1/2 (x,q) = / e^ z p(x,z)dz (28) 

J region 1 /2 

and choose the two regions so that they have no overlap. 
Then 

(pi(x,q)h(0,-q))=2cos{xQ){A lesl (q)Al eg2 (q)) (29) 

where we defined A lcgi (q) = J rogion ; e lqz tp\ (z)tp 2 (z)dz 
and we have 



(A Icgl (q)Al K2 (q))^L 



(30) 

Note that because of the different regions the above quan- 
tity does not have to be real. We define 



which contains the Fourier transform of the function 
d'(r\) = ^ sinh 2 r,,, which is a function peaked around 

t] = 1 with height sinh 2 ( 2L y^) and width 8/2. Fourier 
transforming it will give an oscillating part in the 
correlation function. 

Directly calculating \A q \ in the same way as for bosons 



MA cg2 {q)) 



(31) 



This quanity as well as the correlation function is shown 
in Fig. [5] where the red lines in the plot of the correla- 
tion function in mark Zb og i n and z cn d- | .2 (<?) | 2 shows 
two peaks corresponding to charge and spin velocities. 
The difference in size of the peaks is simply due to the 
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FIG. 4: (Color online) Observable quantity (\A q \ 2 ) and correlations function for parametrically driven two component one- 
dimensional Bose gas. The parameters here chosen are S = .1 and number of oscillations in driving n = (^r) = 20. 



different relative change in the sound velocities u, which 
are given by Eq. ( 10b I , where the relative change in g is of 
course equal for both modes. There are more excitations 
in the spin sector. 



V. CONCLUSIONS 

We have studied interference of parametrically driven 
one-dimensional systems of both bosonic and fermionic 
ultracold systems, and found that the interference pat- 
tern contains information about the excitations driven 
into the systems. For two component systems two kinds 
of excitations are driven on resonance: those correspond- 
ing to charge velocity, and those corresponding to spin 
velocity. This can be extracted from the interference pat- 
tern. In the case of fermions one needs to use different 
regions of the interference pattern to resolve the effect. 
We argued that such experiments could be used to probe 
properties of one-dimensional systems, as measuring the 
sound velocity, and spin charge separation. 
Note added: During the preparation of this manuscript 
the author became aware of Refs. [TH1 HZ], who studied 
similar effects for fermionic systems and obtained similar 
results. 
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Appendix A: Derivation of equation 1 

When we write down the creation operator 
^i. .(xi, zi) for a particle after expansion at posi- 



tion (x\,zi), this particle can have come from either 
quasicondensate, and there is a relative phase between 
those paths 



(Al) 



where ipj(z), j — 1, 2, is the creation operator for a par- 
ticle in quasicondensate j before expansion, and we fac- 
tored out a common phase. The relative phase comes 
from the expansion the particles initially confined in a 
harmonic potential (particles only occupy the ground 
state of the transversal confinement). Consider first a 
particle in a harmonic oscillator at position x ± d/2, 
whose potential is suddenly removed. The wave function 
f±(x, t) expands and picks up a phase factor 



f±(x,t) = 



where i?n = 



1 



T 3 /*y/K 



exp 



(x± 



1 



J mRi 



2R\ 



is the initial radius of the Gaussian 



wave function, m is the particle mass, and wj_ is the har- 
monic oscillator potential, here the transverse confine- 
ment frequency, and 



H, — K n 



M 

mRr 



For large expansion times, so that Rt 3> Rq, Rt 
and 



f±(x,t) 



1 



^I A JR; 



exp 



m ( d 
Itzt- \ x ± — 
2ht \ 2 



(A2) 



tuRq '■ 



(A3) 



As the two systems are initially independent, 
(i/' 1 (z)'!/>2(z)) = 0, and thus the expectation value 

Wioud( x ' ^cioudO, z)) = const. (A4) 




FIG. 5: (Color online) Fermions with spin: (a) (\Ai : 2(q, T)\ 2 ) shows a peak for the spin mode and one for the charge mode. To 
obtain these peaks, one has to integrate over different regions of the interference patter, see text, (b) Correlation function for 
driven two component system of fermions, and integration region. 



However, each single shot will have interference fringes 
due to interference of the different paths, which one can 
see in the expectation value of the density-density cor- 
relation function (p(xi, Z\)p{x2, 22)) of the cloud. The 
operator for finding one particle at r\ — {x\, z\) 1 and the 
other one in ri = (x2, 22), is proportional to 

4ioud(ri)4ioud(r 2 ) oc Mzi)Mz2)e- l{xi+X2)Q/2 

+Mzi)Mz2)e- tAxQ/2 
+Mzi)Mz2V AxQ/2 

(A5) 

here Ax = X\ — X2- The first two terms describe the 
case when the two particles come both from the same 
quasi condensate. The latter two terms describe the case 



when the two particles come each from a different quasi 
condensates: those are the terms giving an oscillating 
term in Eq. (JlJ - and the interference pattern. When 
we multiply the above with its complex conjugate, and 
take the expectation value, the following two terms give 
non-constant contributions: 



(^^{(z^M^M^e 1 ^ + H.c. 
±(tt(z 2 )Mzi))(rt(zi)Mz2)y AxQ + U.c 

±2 $\Z2)${Z{)) 



cos(Air(3) 



(A6) 



the upper (lower) sign is for bosons (fermions). In the 
first step we assumed that systems 1 and 2 are indepen- 
dent, thus the expectation value factorizes, in the second 
step we assumed that they are equal, i.e. in the same 
quantum state, and used ip to describe that state. 
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